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Abstract 

Proton and carbon ion therapy is an emerging technique used for the treatment of 
solid cancers. The monitoring of the dose delivered during such treatments and the 
on-line knowledge of the Bragg peak position is still a matter of research. A possible 
technique exploits the collinear 511 keV photons produced by positrons annihilation from 
P + emitters created by the beam. This paper reports rate measurements of the 511 keV 
photons emitted after the interactions of a 80 MeV/u fully stripped carbon ion beam at 
the Laboratori Nazionali del Sud (LNS) of INFN, with a Poly-methyl methacrylate target. 
The time evolution of the /3 + rate was parametrized and the dominance of 11 C emitters 
over the other species ( 13 A", 15 0, 14 0) was observed, measuring the fraction of carbon 
ions activating /3 + emitters to be (10.3 ±0.7) ■ 10 -3 . The average depth in the PMMA of 
the positron annihilation from /3 + emitters was also measured, Da+ = 5.3 ± 1.1 mm, to 
be compared to the expected Bragg peak depth Dsragg = H-0 ± 0.5 mm obtained from 
simulations. 
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Introduction 

The use of proton and carbon ion beams has become more and more widespread as 
an effective therapy for the treatment of solid cancers (hadrontherapy) . These beams 
have maximum energy density released at the Bragg peak (Fig. [I] Top) at the end of 
their range, as opposed to X-rays or 7-rays, which are absorbed by the body with an 
exponential decrease in the delivered dose with increasing tissue depth after an initial 
build-up pp. 

Due to their very favorable profile of released dose in tissue, the charged hadron 
beams can be very effective in destroying the tumor and sparing the adjacent healthy 
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Figure 1: (a) Dose distributions as functions of water-equivalent depth estimated with the 
GEANT MC for 279.2 MeV 12 C nuclei in water (solid-line histogram), PMMA (dashed his- 
togram) and bone (dotted histogram). Experimental data for depth dose distributions in PMMA 
are shown by points pQ. (b) Distributions of positron-emitting nuclei produced in these materials 
as a function of water-equivalent depth [T]. 



tissue in comparison to the standard X-ray based treatment [2]. On the other hand, the 
higher spatial selectivity of hadrontherapy asks for a dedicated approach to the delivered 
dose monitoring. 

The uncertainty on the position of the dose release in hadrontherapy treatment can 
be due to different factors, i.e. quality and calibration of the Computed Tomography 
(CT) images, possible morphologic changes occurring between CT and treatment, patient 
mis-positioning and organ motion during the treatment itself. All these effects give an 
overall uncertainty of the order of few millimeters, that can be larger than the dimension 
of the dose release spot at the Bragg peak. 

Several methods have been developed to determine the Bragg peak position by ex- 
ploiting the secondary particle production induced by the charged hadron beam, and 
among these one of the most promising is the PET-like technique: the collinear 511 keV 
photons produced by positrons annihilation from j3 + emitters created by the beam are 
measured. The relationship between the (3 + emitters densities and the dose release has 
been studied with Monte Carlo (MC) simulations, as shown in Fig. [I] [l a . The measure- 
ment of the rates of such emitters can also provide precise monitoring of the dose, which 
is in turn essential for a good quality control of the treatment. This technique has already 
been developed with measurements both after the irradiation [31 2] and on-beam [31 [5] • 
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To this aim, measurements are needed to allow the MC tuning, which is critical for the 
appropriate development of the technique. In particular, a determination of the isotopic 
composition of the f3 + emitters and the corresponding rates has not been performed 
for the carbon ion treatments. Furthermore, papers in literature mostly report time 
integrated measurements and do not investigate the time structure of the emission in 
presence of irradiation. 

In this paper we present measurements of the properties of the /3 + emitters, produced 
in a Poly-Methyl Methacrylate (PMMA) phantom during an irradiation with carbon ions, 
by observing the two 511 keV photons produced in the positron annihilation (7— PET). 
In section [lj we describe the setup of the experiment performed at the INFN Laboratori 
Nazionali del Sud (LNS) in Catania in the interaction of a 80 MeV/u fully stripped carbon 
beam with a PMMA target (see [7] for a related experiment with the same setup) . Data 
analysis tools are detailed in Sec. [2] 

With the collected data we investigated three aspects of the /3 + emission: the isotopic 
composition of the emitters (Sec. |3.1[ ), their corresponding rates, via a study of the time 
dependence of the emissions during irradiation (Sec. 3.2), and the position of the /3 + 
emitters with respect to the Bragg peak (Sec. [4]). With respect to the existing studies 
we have allowed for incident beam intensity variations. 




Figure 2: Schematic view of the experimental setup: the Nal detectors are placed 20 cm away 
from the PMMA. The acquisition is triggered by the coincidence of the two crystals. On the 
right the picture of the actual experimental setup is shown. 



The experimental setup is shown in Fig. [2J A 4 x 4 x 4 cm 3 PMMA target is placed 
on an 80 MeV/u fully stripped 12 C ion beam. The beam rate, ranging from hundred 
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of kHz to ~ 2MHz, is monitored with a 1.1 mm thick scintillator on the beam line, 
placed at 17 cm from the PMMA, and read-out by two photomultiplier tubes (PMTs) 
(Hamamatsu 10583) put in coincidence (Start Counter). 

A pair of cylindrical sodium iodide activated with thallium crystals Nal(Tl) (r = 
2.5 cm and h = 5 cm) is placed at 45° (225°) with respect to the beam line, at 20 cm 
from the PMMA. The scintillation light of the two crystals is detected by two Scionix 
VIA— EI PMTs triggered in coincidence within a time window of 80 ns. A 12-bit QDC 
(Caen V792N) and a 19-bit TDC (Caen V1190B) provide the measurements of energy 
and arrival time. We use Nal crystals to detect the 7— PET signals because of their 
high light yield and energy resolution in the 0(1 MeV) ranges. In order to perform 
position-dependent measurements, the PMMA can be moved along the beam axis with 
an accuracy of 0.5 mm. 



2. Data selection 



From the charge collected in the two Nal detectors we measure the energy of the two 
photons (E*, i = 1, 2) from the positron annihilation. To calibrate the Nal detectors we 
used 22 Na (511 keV), 137 Cs (662 keV) and 60 Co (1.17 MeV and 1.33 MeV) sources. A 
good linearity is obtained up to 1.5 MeV. The relative energy resolution at 511 keV is 
a E /E = 2.7%; we define the signal window 0.48 < E. t < 0.53 MeV. 

Fig. 3] shows the correlation between the E^ and E 2 and the energy spectrum of 
both detectors after requiring the other one to be in the signal window. A /3 + decay is 
identified by the presence of two back-to-back photons with an energy compatible with 
E = 511 keV. Thus the events with both energy depositions in the signal window provide 
a background free measurement of the number of /3 + emitters. 




Figure 3: Left: correlation between the calibrated energy depositions in the two Nal detectors. 
Right: calibrated energy spectrum in both Nal detectors after requiring one of the two to be in 
the signal region (mean = (0.507 ± 0.012) MeV). 



To evaluate the detection efficiency (tdet) we simulated the experimental setup with 
GATE [8 , a simulation tool dedicated to medical imaging, radiontherapy and hadron- 
therapy based on the GEANT4 MC code 0(10]. 
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In order to measure the position of the peak of the f5 + decay spatial distribution it 
is important to evaluate the dependence of e^ e t on the position of the two photon vertex 
with respect to the Nal detectors. To this aim we define x as the distance along the 
beam line between the two photon vertex in the PMMA and the interception between 
the beam line and the line connecting the two detectors (see Fig. E| . 




D °p + 



Figure 4: Definition of the x and the D variables (see text). 



The resulting functional dependence tdet( x ) is shown in Fig. [5] and was parametrized 
with a single gaussian. 

In order to measure the 7— PET rate, the number of carbon ions Nc must be evalu- 
ated. In a given time interval, Nq is computed by counting the number of signals given 
by the Start Counter (Nsc) within randomly-triggered time-windows of T w — 2 [is. 
From the number of time windows considered (N w ) and the total acquisition time (T to t), 
the number of carbon ions Nc is estimated to be: 

Nc^^r- (1) 
f-sc N w l w 

The Start Counter efficiency esc — (96 ± 1)% has been estimated by exploiting the 
two-sided PMT readout with negligible dark counts. 

The number of impinging carbon ions has to be corrected for the dead-time ineffi- 
ciency, 6dt has been estimated from the total acquisition dead time (Tdead) as: 

-1 Tdead / \ 
(DT = 1 - ■ (2) 



0.0012 




x [mm] 

Figure 5: Detection efficiency as a function of the PMMA position on the beam line {\ 2 /ndf = 
74/24, x mea „ = (-0.046 ± 0.13) mm and a = (13.45 ± 0.09) mm). See text for details. 



The measured values of tor range from 70% at an average carbon ion rate of 0.6 MHz 
up to 47% at 2 MHz. This efficiency correction was then applied to data as described 
below. 

3. Measurements of the time evolution of the /3+ emission 

The rate of (3 + decays and the isotopic composition of the emitters was measured 
as a function of time both during irradiation and in the intervals in between. The time 
dependence of the emission during the irradiation results from two main contributions: (i) 
the creation of new emitters induced by the passage of the carbon ions in the PMMA, and 
(ii) the decay of the previously created ones. When the irradiation time is comparable 
to the decay time of the emitters, the relation between the emitter and dose rates is 
non-trivial. This is the case studied in this paper. 

3.1. Isotopic composition of emitters 

The isotopes that can be produced during the carbon ion irradiation of the PMMA 
are n C, 13 N, ls O, 14 and 10 C with half-lives of 29 min, 15 min, 2 min, 100 s and 19 s 
respectively. The wide spread of lifetimes allows to discriminate among isotopes. We 
chose those isotopes because they are the one with longer half-lives and, as it will be 
seen later on in the "during irradiation" measurements, the gamma distribution is time 
has a rise of the order of tens of minutes. When the beam is turned off (t=0), the time 
evolution of the /3 + decay rate (R(t)), is: 

N B 

R{t)=R(<S)Y j f i e- t ^. (3) 

i=l 
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assuming N s different species of isotopes with fractions /$ and lifetimes 7$. The mea- 
sured 7— PET, expected to be proportional to R(t), in a run with no-beam (measurements 
after irradiation) as a function of time is shown in Fig. p\ 




Figure 6: Distribution of the number of detected 7— PET as a function of time in a run with 
no-beam; the fit with an exponential function is shown (r = 33 ± 3 min, \ 2 /ndof = 15.3/13). 



Since the time t is measured with respect to the start of the run and not to the end 
of the irradiation as required by Eq. [3] the sensitivity to the fast decaying isotopes ( ls O 
and 14 0) and the combination of different runs are compromised. We have therefore 
decided to analyze only the no-beam run with the largest statistics (shown in Fig. [6|. 
The available data does not allow us to perform a fit with a double exponential function 
accounting for both 11 C and 13 iV populations and therefore we can only test the hypoth- 
esis that one isotope dominates. The fit with a single exponential function results in a 
lifetime r = 33 ± 3 min (\ 2 /ndof = 15.3/13), thus indicating that the 11 C component 
is dominant. 



3. 2. Time evolution of the emission rate during irradiation 

As stated in the introduction of section [3J the time evolution of the j3 + emitters 
becomes non trivial when irradiation lasts for a time interval comparable with the lifetime 
of the emitters. We therefore elaborated a model to describe the simultaneous occurrence 
of new activations and decays of the (3 + emitters and used it to fit to the time evolution 
of the measured rates, testing its goodness and extracting its parameters. 
Defining the number of /3 + emitters for each of N s species as iVJ ec (s — l,...,N s ), with 
t s lifetime, we can write: 

d ^=A s d J^-lN !ec (t). (4) 
dt dt t s decK ' K ' 

The first term on the right side of the Eq. [4] represents the activation induced by 
the impinging carbon ions, with A r s aw representing the fraction of carbon ions activating 
/3 + emitters. The second term takes into account the decay of the activated nuclei with 
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decay times r s . Each species has an independent time evolution, so the measured 7- 
evolution in the detector (iV 7 ) takes into account all species contributions: 



dN 7 (t) 
dt 



N N s 

tdet > — 

s — 1 



(*) 



PET 



(5) 



Eq. H requires the knowledge of Nc(t). As it can be seen for two different runs shown 
in Figs.Tf](Left) and [9] (Left), Nc{t) behavior cannot be always described with a function 
that allows Eqs. [1] to be solved. We therefore developed two different methods: (i) an 
analytical one to be used when Nc(t) can be described by an analytical function, and 
(ii) a numerical one that can always be applied, but requires the assumption of one 
dominant isotope among the emitters. It is to be noted that accounting for the time 
dependence of the beam intensity is particularly needed when the irradiation lasts for 
periods comparable to the decay times of the isotopes, as it is the case in the therapeutical 
treatments. 





20 30 40 50 



Figure 7: Measured N c (Left) and N 1 (Right) as a function of time. The fit of the measured JV 7 
distribution using Eq. [7] (solid line) gives the species ratio 11 C over 13 N (see Eq. [§]). 



3.2.1. Analytical method 

The Fig. [7] (Left) shows the measured Nc(t) fit by an exponential function (Nc(t) = 
Vbeam exp(— t / Tb eam )) . In this case, the solution of Eqs. [JJand [5] is: 

Nj = (det y ] I I 1 \ I "beam ^ + — e~ ^ . (6) 




TV, 



with 



-£det v be 



A s 

1 



U l)( 



In Sec. [XT] the lifetime measurement obtained from the fit of Fig. [6] indicates that the 
C component dominates. By using Eq.§ we can now measure the relative contribution 
to Nry from 11 C and 13 iV species, Ac and Am- Assuming to as the time for which 



Ng c (t ) = N* c (t ) = 0, Eq. [6] can be written as follow: 



8 




iV 7 = e det 74jv — J h J 1 Vbeam e T <>">™ - e T 6eam . (7) 



1 fair 



Fitting the time distribution of N 7 with Fig. [7] (Right), we obtain: 

—^ = 16.6 ±2.7. (8) 

A N 

The dominance of 11 C over 13 N is therefore confirmed and the measured ratio is 
consistent with the results, obtained with a higher energy beam, in Ref. [5J. 

3.2.2. Numerical method 

In general the carbon ion rate cannot be parametrized with a simple function allowing 
an analytical solution of Eq. [4] Therefore, a numerical method valid in the hypothesis of 
one dominant isotope only has been developed. Defining the integral of the number of 
carbon ions (N c ) and of the detected 7— PET (Nl) in the time bin i with E>w width, we 
can estimate, by approximating the derivative with the differential increment, the A r 
parameter in each single bin: 

where the raw suffix indicates that the measurement has not been corrected for the 
detector efficiency. As an example, the numerical method has been applied to the Nc 
and Nj measurements shown in Fig. [7] The measured A r i aw parameters using Eq. [9] are 
reported in Fig. [8j By minimizing the \ 2 > defined taking into account correlations and 
assuming that the measured parameter is independent of the time bin, we obtain a mean 
value A™ w = (4.7 ± 1.0) • 10~ 6 , with x 2 /DOF = 37/69. 

Conversely, we can exploit the same equations Q, assuming a measurement of A™ w 
from a calibration run, to estimate the number of incident carbon ions (N^,), once the 
numbers TV* (Eq. |9j) are measured for a specific data sample: 

C A rau, y e * DT Bw y e l jjrp ' ' ' 1 





We have tested this procedure on a long run (see Fig. [9] for the measured N 7 and 
Nc) with slightly unstable beam conditions. The A™ w parameter was estimated on the 
shorter run from which Fig. [8] was taken, after adjusting for the detector efficiency as 
detailed in the next section. 



Fig. 10 shows the cumulative distribution of the number of ions estimated by the 
measurements of 7 — PET, N c , blue dot points, compared to the cumulative number 
of carbon ions measured with the Start Counter N^, c , magenta dashed data. A good 
agreement is visible, also at times comparable with the lifetimes of the decaying isotopes. 
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Figure 8: A raw parameter time evolution and fit to a constant. 




Figure 9: Measured N c (Left) and 7V 7 (Right) as a function of time. 



4. Spatial dependence of the /3+ emission 

The proposed numerical method to estimate Nc relies on the knowledge of the A r n 



parameter (see Eq. 10 ). The A™ w parameter measurement described so far was corrected 
only for the dead-time efficiency. To obtain the final A parameter we must correct for 
the detector efficiency {edet)- The latter depends on the position of the f3 + emitters along 
the beam line with respect to the line connecting the two Nal detectors, as discussed in 
Sec. [2] In order to study the spatial distribution of the emitters the A™ w parameter was 
measured for several positions of the PMMA, i.e. for several values of D as defined in 
Fig. [4] It can be written as: 

/r] A 
e d et(x) -£-(x\D)dx (11) 
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Figure 10: Cumulative distribution of the number of carbon ion measured with the Start Counter, 
N^ c (magenta dashed data), compared to the number of ions estimated by the measurements 
of 7 — PET, (blue dot points). The plot on the right is a zoom in the first 50 min of 
acquisition. 



with dA ec /dx(x\D) the density of emitters as a function of the depth x in the PMMA 
(as defined in Fig. [4]) . The dependence on D takes into account the distribution shown 
in Fig. [I] and the relative position of the PMMA and the Nal detectors. Since the litera- 
ture [T] reports widths of the /?+ emitter peak distribution smaller than e(x) (comparison 
of Fig. [ljand Fig. [5]), we approximate dA ec {x\D) / dx ~ 5[x — (D — Dp+)) where Dp+ is 
the position of the j3 + emitters within the PMMA. Finally we obtain: 

AZ w {D)=A Q -e det {D-Df j+ ) (12) 

with A the constant component of the A parameter, D the distance along the beam line 
between the line connecting the Nal crystals and the beam entrance face of the PMMA 
and Dp+ the mean position of the emitters in the same reference system, 
shows the dependence of the measured A parameter from D: the width of 
the distribution is consistent with the one of the efficiency in Fig. [5] thus confirming the 
hypothesis used to obtain Eq. p~2| 

From the fit of the spectrum in Fig. [TT] with the function in Eq. [12] we obtain: 

A = (10.3 ± 0.7) -10" 3 (13) 
D p+ = (5.3 ±1.1) mm. (14) 



(see FigJI 
Fig. In 



The dose deposition in the PMMA has been simulated with FLUKA [TT] . The result 
is shown in Fig. 12 as a function of the depth from the beam entrance face of the PMMA 
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Figure 11: The measured A™" (D) parameter as a function of the position of the PMMA, fitted 
as described in the text. 



with the beam entering form the left side of the plot. In our configuration the Bragg 
peak is at 11.0 ± 0.5 mm from the beam entrance face (light yellow band) of the target. 
The simulation is confirmed by the picture of the PMMA after the data taking , shown 



in Fig. 12 (Inset): the dose distribution is visible by the deterioration (light yellow band) 



of the PMMA. 




Figure 12: Dose deposition in the PMMA simulated with FLUKA as a function of the depth 
from the beam entrance face of the PMMA. 
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Using the measured value of Dp + (Eq. 14) with the simulated Bragg peak position 
D Bragg j we obtain the distance £ between the /3+ emitters and the Bragg peaks: £ = 
|%-%, M | = (5.7±1.2)min. 



Conclusions 

We presented a study of the rate of 7— PET produced in the interaction of 80 MeV/u 
fully stripped carbon ions with a PMMA phantom, using a pair of Nal crystal detectors. 

We proposed, and validated with data, a model to describe the activated nuclei j3 + 
decay during the irradiation. We demonstrated the possibility to estimate the number 
of impinging carbon ions from the number of observed 7— PET. 

We measured the ratio between the number of activated 11 C and 13 N to be Ac /Am = 
16.6 ± 2.7 and a number of (10.3 ± 0.7) • 10~ 3 generated 11 C ions, per impinging carbon 
ion undergoing /3 + decay. 

Finally we measured the mean position of the /3 + emission to be Dp + — (5.3±1.1) mm 
from the beam entrance to the PMMA face, to be compared to the simulated Bragg peak 
position DBragg = (11.0 ±0.5) mm. Of course more accuracy on the measurement could 
be achieved with larger scale experiments, but anyhow we think that fair information 
can also be achieved from the comparison with the presented data. 

All this information can be used as a benchmark for the f3 + emitters MC simulation 
of hadrontherapy. 
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